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Motivated by recent synthesis of the hyper-honeycomb material y 5 -Li 2 lr 03 , we study the dynamical struc¬ 
ture factor (DSF) of the corresponding 3D Kitaev quantum spin-liquid (QSL), whose fractionalised degrees of 
freedom are Majorana fermions and emergent flux-loops. Properties of this 3D model are known to differ in 
important ways from those of its 2D counterpart - it has flnite-temperature phase transition, as well as distinct 
features in Raman response. We show, however, that the qualitative behaviour of the DSF is broadly dimension- 
independent. Characteristics of the 3D DSF include a response gap even in the gapless QSL phase and an 
energy dependence deriving from the Majorana fermion density of states. Since the majority of the response 
is from states containing a single Majorana excitation, our results suggest inelastic neutron scattering as the 
spectroscopy of choice to illuminate the physics of Majorana fermions in Kitaev QSLs. 


Introduction. A central theme in an ongoing search for new 
states of matter is the interplay of interactions and dimension¬ 
ality, leading to new exotic phases displaying e.g. topologi¬ 
cal order and quantum number fractionalisation. Prime ex¬ 
amples are magnets in which quantum fluctuations (QF) are 
strong enough to suppress long-range magnetic order, giving 
rise instead to quantum spin-liquids. One route to these is 
via reduced dimensionality, as shown by many experimen¬ 
tal examples of QSLs in spin chain systems. For such ID 
systems, the materials realisation and theoretical understand¬ 
ing of spin-liquid behaviour is well established [1, 2]. In 
higher dimensions an alternative route is required, and An¬ 
derson in his original proposal [3] of resonating valence bond 
(RVB) states recognised the importance of frustration. Here, 
however, physical examples are scarce, and theoretical under¬ 
standing is less complete. Indeed, it took several decades be¬ 
fore even the existence of RVB phases for microscopic 2D 
model Hamiltonians was firmly established theoretically [4] . 

Theoretical tools and solvable models are important for de¬ 
veloping the understanding of unconventional phases. A num¬ 
ber of powerful methods exist for ID, such as density matrix 
renormalisation group (DMRG) [5] and Bethe-Ansatz, which 
allow for a quantitative comparison between theory and ex¬ 
periment [6]. For 2D and 3D QSLs, a qualitative understand¬ 
ing has been developed using large-W limits [7, 8], appro¬ 
priate mean held theories [9] and the identification of long- 
wavelength descriptions [10, 11]. These approaches do not, 
however, provide detailed results for specific realistic models. 

Against this background, the Kitaev honeycomb model [12] 
occupies an important position as a rare instance of an exactly 
solvable 2D model that supports a variety of different QSL 
phases, stabilised by exchange frustration. Moreover, its 3D 
generalisations [13] provide an opportunity to examine the in¬ 
fluence of dimensionality on these phases. 

Materials that are candidates for the realisation of Kitaev 
frustration should have dominant spin-orbit coupling, since 
exchange interactions in the model link together real-space 
and spin-space anisotropies. Recently, a set of 2D honeycomb 
lattice compounds, {Na,Li}2lr03 iridates [14-18] and RuCls 


[19-24] have been proposed to exhibit dominant Kitaev-like 
spin exchange arising from spin-orbit interactions [25]. 

In addition, in two insulating 3D polymorphs (3- and 7- 
Li2lr03 [17, 18], the Ir^+ have been shown to form struc¬ 
tures dubbed hyper-honeycomb and stripy-honeycomb lat¬ 
tices. These are part of a whole series of “harmonic" hon¬ 
eycomb lattices which might realise Kitaev QSL physics 
in 3D via the interplay of spin orbit coupling, interactions 
and the simple fact that each lattice site has coordination 
number three [17, 26-30]. Recent work [31] has identified 
Jeff = 1/2 degrees of freedom in the low-energy description of 
;5-Li2lr03 and furthermore confirmed that the effective spin- 
model has dominant Kitaev-like exchange. So far, all experi¬ 
mental candidate materials show long-range magnetic order at 
low temperatures. However, an observation of non-coplanar 
spiral magnetism in 3D polymorphs is in itself a strong in¬ 
dication of dominant Kitaev exchange [32-34], and perhaps 
a proximate QSL. Also, it has been shown that the QSL of 
the Kitaev model is stable with respect to small integrability- 
breaking Heisenberg perturbations, while larger perturbations 
result in various magnetically ordered states [26, 27, 31]. 

Here, we present the results of an exact calculation of the 
dynamical spin response for a 3D QSL, the Kitaev QSL of the 
hyper-honeycomb lattice [13] of /3-Li2lr03. Our results show 
signatures of spin fractionalisation into emergent quasiparti¬ 
cles, Majorana fermions and flux-loops. The character of the 
QSL phase itself reflects the dimensionality of the lattice. In 
particular, point like fluxes of the honeycomb lattice in 2D are 
replaced by extended flux-loops in 3D, which leads to a true 
flnite-temperature phase transition [35]. 

Our central finding is that the dynamical structure factor 
(DSF) is a rather direct probe of the Majorana fermion excita¬ 
tions in Kitaev models, independently of dimension. While 
details vary - in 3D, fine structure from several Majorana 
bands is present, while in 2D features that stem from van 
Hove singularities in the Majorana density of states (DOS) 
are more prominent - the changes are small. This is in con¬ 
trast to other dynamical probes, such as Raman scattering, 
which shows qualitative differences between 2D and 3D Ki- 
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taev QSLs, e.g. via polarization dependence [36, 37]. 

The DSF measures the excitations induced in the ground 
state by a spin flip. In Kitaev models these consist of static 
fluxes and a variable number of Majorana excitations. This 
perspective has already enabled us to develop a complete theo¬ 
retical picture for a dynamical response in all different gapped 
and gapless, Abelian and non-Abelian QSL phases of the 2D 
honeycomb Kitaev model [38, 39]; in fact, signatures of the 
Majoranas have arguably already been observed in the short 
time dynamics (high energy response) of RuCls compounds 
[22, 24]. Calculations for the 3D model present a signiflcant 
additional challenge, made feasible only by the technical de¬ 
velopments we have described elsewhere [39]. 

Model. The hyper-honeycomb 3D Kitaev model is built 
from spin- 1/2 degrees of freedom arranged on a hyper¬ 
honeycomb lattice (see Fig. 1) that interact via bond- 
dependent, nearest-neighbour Ising exchange Ja. Three lat¬ 
tice directions are labelled by a = x^y^z referring to compo¬ 
nents of spins involved in the exchange bond a. In terms of 
Pauli matrices and using the notation {jk)a to indicate two 
sites, j and k, connected via a bond a, the Hamiltonian reads 

^ = -y - y - Y ^j^k- (i) 

{jk) X {jk)y {jk) Z 

The model has two types of ground state, namely gapless and 
gapped QSLs for | Jq, | ^ | J/ 3 1 + | /y | (with a ^ [3 ^ ^ ^ a) 
[13], and the topology of the phase diagram is the same as 
in a number of other 2D and 3D Kitaev models [12, 26, 29]. 
We note that the Majorana spectrum in 2D honeycomb/3D 
hyper-honeycomb possesses a single Dirac point, and a gap¬ 
less nodal line respectively. 

Our calculation has as its starting point the re-expression of 
this Hamiltonian, following the original approach of Kitaev 
[12], in terms of Majorana fermions. Of this we briefly men¬ 
tion the features important for understanding the physical con¬ 
tent of our results; we discuss what is technically new in three 
dimensions below, and refer the reader to the methods we have 
developed in the context of two dimensions Refs. [38, 39]. 



Figure 1. Schematic picture of the unit cell of a hyper-honeycomb 
lattice showing three inequivalent bond types x,y,z together with 
lattice vectors ay Yellow sites make up the primitive unit cell con¬ 
sisting of four Ir atoms in / 3 -Li 2 lr 03 . 


The spin on site j is expressed in terms of four Majorana 
fermions: Cj^h^, with a = x^y^z. The spin operators can 
then be represented as = icjh^. The two ‘species’ of Ma¬ 
jorana fermions take on separate roles. The 6 ’s combine to 
yield emergent non-dymanic fluxes. These are captured via 
loop operators W/; see Refs. [13, 28] for their precise def¬ 
inition in terms of spin operators. The loop operators have 
eigenvalues W/ = ±1 associated with the Z 2 -fluxes pierc¬ 
ing the irreducible loops. These fluxes are point-like excita¬ 
tions in 2D, but form lines in 3D. The b Majorana fermions are 
complemented by the c’s, which can be combined in pairs to 
make standard complex dynamical fermions. Using the con¬ 
served quantities corresponding to loop degrees of freedom, 
the Hilbert space can be decomposed into a ‘flux’, \F) and 
a ‘matter’, |M), sector. For a given flux sector, a Hamil¬ 
tonian for the matter fermions describes hopping between 
sites in the presence of the Z 2 gauge fluxes. In the follow¬ 
ing, we denote the ground state of H (which is flux-free) as 

|0) = |Fo)^|Mo). 

Dynamic structure factor. Our central result is the numer¬ 
ically exact evaluation of the DSF 5'(q, uj), deflned as 

= (*)> ( 2 ) 

ajk 

where = ( 0 |(Jj (t)(j^( 0 )| 0 ) is the time-dependent spin 
correlation function. The DSF is directly related to the cross- 
section measured in inelastic neutron scattering (INS) exper¬ 
iments [40], and at q = 0 to the signal obtained in ESR ex¬ 
periments. At the heart of our analysis is the observation that 
the static flux degrees of freedom (the 6 ’s) can be eliminated 
from the calculation which leads to strictly zero spin correla¬ 
tions beyond nearest neighbours and a local quantum quench 
of Majorana fermions [41]: a Majorana fermion is added to 
the ground state \Mq) of Hq which then evolves under the ac¬ 
tion of a matter Hamiltonian corresponding to a different flux 
sector. For each component this flux sector is obtained 
by flipping the Z 2 gauge held on the corresponding a-lattice 
bond (ij)a- Note that the increased number of inequivalent 
bonds in the hyper-honeycomb lattice compared to the hon¬ 
eycomb one entails an increased number of inequivalent local 
quench problems entering Eq. 2. Despite the exact solvabil¬ 
ity of the model, the non-equilibrium nature of the response 
makes it notoriously difficult to calculate the DSF. While in 
2D one can resort to flnite-size numerical calculations [39], 
these become useless for 3D. Using the exact semi-analytical 
approach [38] which we developed in the studies of the 2D 
case we have been able to solve the problem also in 3D in the 
thermodynamic limit. 

The Lehmann representation of the DSF provides a number 
of qualitative insights. Let {|A)} be the many-particle eigen¬ 
states of the matter Hamiltonian in the presence of four flux- 
loops with a common a-bond. Then 

= 2nFy ^(Mo|c,-|A)(A|c,|Mo)(5(a; - Ay), (3) 
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Figure 2. Dynamical structure factor as a function of ujI Jz along four 
high symmetry cuts in the Brillouin zone, shown for the isotropic 
point Jx = Jy = Jz on a linear colour scale. Note that in order 
to improve the visibility of the lower intensity features especially at 
high energies all intensity values above 30 (arb. units) are shown in 
dark red. For the full range of intensity see Fig. 3. 


where G — i} depend on the spin component, and 
is the excitation energy, which includes the contribution 
from the flux-gap. This representation reveals the main quali¬ 
tative features of the response. They are remarkably indepen¬ 
dent of the dimensionality. First, the response vanishes below 
the flux-gap, which is the energy cost of exciting flux-loops by 
flipping a bond. Strikingly, this response gap is present even 
in the QSL phase with gapless Majorana excitations. Note 
that the value of the flux-gap depends signiflcantly on the spin 
component probed if the J^’s are not all equal. The gap is the 
most direct consequence of fractionalisation, and is the signa¬ 
ture feature of the emergent Z 2 gauge fleld. The response 
above the flux-gap is entirely due to matter fermions. Al¬ 
though it includes multi-particle contributions, we And that the 
single-particle sector is dominant (accounting e.g. for 87% of 
the response at q = 0 with all Ja equal). INS therefore probes 
the Majorana fermion DOS rather directly. 

Results. It was suggested recently [31] that the effective 
Hamiltonian for ; 5 -Li 2 lr 03 is likely to be in the regime close 
to the ferromagnetic isotropic point, and so we concentrate 
on this point in parameter space in the following. We note 
that behaviour is different in the opposite limit of one domi¬ 
nant exchange. For example, the component of response 
corresponding to the dominant interaction Ja has a sharp 5- 
function contribution from a pure flux excitation whose pres¬ 
ence/absence gives rise to a dynamical phase diagram similar 
to the 2D case [38] (not shown). 

Fig. 2 shows the cc-dependence of the DFS along cuts in the 
Brillouin zone (BZ). The response vanishes below the flux- 
gap, although this is a gapless QSL phase. Gross features of 
DSF are remarkably similar to those of the 2D honeycomb 
case [38]. This is due to a similar linearly vanishing Majorana 
DOS both on the 2D honeycomb and 3D hyper-honeycomb 


lattice, and the fact that the DSF originates from a local quan¬ 
tum quench despite the extended nature of the flux-loops in 
3D. Response above the flux-gap varies smoothly with energy, 
and falls off rapidly above the fermion bandwidth. Extra flne- 
structure due to the 3D nature of the lattice is generated from 
multiple Majorana bands originating from the increased size 
of the unit cell. The momentum dependence of the DSF in 
the Brillouin zone for three flxed values of ujj Jz is shown in 
Fig. 4. The limited variation reflects the fact that spin correla¬ 
tions in the QSL are short-range. We And that there is a strik¬ 
ing “inflexion point” at energy ujj Jz ~ 2.2 across which the 
relative intensities in the BZ reverse. At the inflexion point, 
the response is essentially independent of wavevector across 
the BZ. The reversal can be traced back to mixing of differ¬ 
ent contributions from nearest-neighbour (n.n.) and same-site 
(s.s.) correlators as shown in Fig. 3. Note the characteristic 
peak (dip) for the s.s. (n.n.) components around uj = 2.1 Jz, 
which is related to the corresponding peak in the DOS as 
shown in the inset. Furthermore, at the edge of the BZ there 
is dependence only on the Xi — T, T — Ai direction which 
is true for all values of uj. These features are also independent 
of small anisotropies of the Ja. 

Finally, the results shown in Fig. 5 compare the powder av¬ 
eraged DSF for the ferromagnetic and the antiferromagnetic 
isotropic points. The notable distinction is that the intensity 
peak shifts from |q| = 0 for the former to the boundary of 
the BZ (around |q| « 3) for the latter. Compared to the hon¬ 
eycomb lattice case [24] we And three, instead of two, broad 
modes in energy. 

Discussion. The presence of a flux-gap and the broad re¬ 
sponse at higher energies provide measurable signatures of 
fractionalized spin-excitations - emergent gauge fluxes and 
Majorana fermions - in INS experiments. While / 3 -Li 2 lr 03 
materials are long-range ordered at low temperatures, it might 
be nevertheless possible to observe the signatures of fraction¬ 
alisation in dynamical scattering experiments. This should 
hold as long as emergent quasiparticles are only weakly con- 



Figure 3. The four inequivalent components of the dynamical struc¬ 
ture factor at the isotropic point Jx — Jy — Jz- Inset: band resolved 
DOS of the Majorana fermion dispersions. 
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Figure 4. Dynamic structure factor for Ja = 1, a = x,y,z ai fixed values of (jjjJz. Intensity is shown for high symmetry planes in the 
Brillouin zone. The left and right figures show complementary momentum dependence while the center figure shows an ‘inflexion point’ at 
the transition between the two. Left and right figures are representative of the response below/above the “inflection point” aiuj/ Jz ~ 2.2 


fined, so that they represent the natural degrees of freedom to 
account for the short-time behaviour, which is governed by a 
dominant Kitaev interaction. 

Despite fundamental differences between the Kitaev QSL 
phases in 2D and 3D (point like fluxes versus loops, cross-over 
versus true phase transition at finite temperature) we find that 
increased dimensionality does not affect the qualitative be¬ 
haviour of the dynamical structure factor. The latter is surpris¬ 
ingly similar on the 2D honeycomb and 3D hyper-honeycomb 
lattices, in contrast to other response functions such as Ra¬ 
man scattering [36, 37], which shows qualitative differences 
between 2D and 3D. This observation suggests that INS is 
the method of choice for a spectroscopic probe of Majorana 
fermions. The way INS couples to these is sufficiently com¬ 
plex to provide a picture of their full spectrum, but simple 
enough to allow for a direct measure of their properties. 

Indeed, the insensitivity of spin-correlations to dimension¬ 
ality suggests a great level of universality for Kitaev QSLs on 
the harmonic honeycomb series. This is related to the fact 
that on all these polymorphs the dynamical response origi¬ 
nates from a local quantum quench, which arises from the hi¬ 
erarchy of the fractionalized quasiparticles - gapless Majorana 
fermions and static fluxes. 

This viewpoint in terms of emergent degrees of freedom 
and their interplay in a local quench setting [41] thus appears 
as the natural language for probing fractionalisation in Kitaev 
spin liquids. It immediately suggests that instead of changing 
dimensionality, qualitative differences will rather result from 
a change of the low energy DOS. In particular, studying the 
DSF for a QSL with entire gapless zero energy surfaces [29] 
or only Weyl points [30] is an obvious and interesting subject 
for future research, which can be tackled with the methods 
we have developed. More ambitiously, it will be of interest to 
see how magnetic instabilities - in particular those resulting 
from integrability breaking interactions - will dependent on 
dimensionality on one hand, and the properties of the matter 
fermions on the other. Such studies will be necessary to make 
detailed, fully quantitative contact with experiment. How¬ 
ever, in the presence of integrability-breaking perturbations, 
the methodological situation - in the absence of the exact ap¬ 


proach we have developed for the pure Kitaev models - seems 
to be considerably more daunting. 
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